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INTRODUCTION 


Thermal convection is a phenomenon which is very common in /*3 
astro- and geophysical conditions. Energy transport by convect- 
ion occurs in the nuclei of stars of the upper part of the main 
sequence and in the shells of the stars of Its lower section, 
particularly on the Sun. The interaction of convection with 
radical pulsations and rotations is apparently the mechanism 
which causes the cepheid variable. Let us assume that the band 
structure of the surface of the planets, giants Jupiter and Saturn, 
is the reflection of the convective instability of the atmospheres 
of these planets. On the Earth, convection influences the global 
motion in the atmosphere and the oceans, and convection in the 
mantle of the Earth is responsible for the motion of the continents. 


* Numbers in the margin indicate pagination of original foreign text. 
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I 

Many studies have investigated convection under astro- 
and geophysical conditions. However, the problem is so complex 
that we can only think about solving certain particular problems. 

The classical problem of Raleigh-Benar has been known for 
a long time. It covers the convective stability of an infinite 
horizontal layer of an incompressible liquid, in which a vertical 
temperature gradient is maintained 111. However, the convection 
occurring under the conditions in which we are interested is a /4 

much more complex phenomenon. It is basically related to the fact 
that convection is not an individual, isolated process, but occurs 
in interaction 5 with other phenomena of stellar or planetary dynamics 
and thermodynamics. Thus, in many astrophysical cases, it is neces- 
sary to consider the influence of rotation and magnetic fields upon 
convection. Layers which are convectively unstable in real objects 
usually are connected with thermally stable layers. The penetra- 
tion of convective pulsations in these layers must be considered 
in order to obtain the correct physical picture. 

Frequently, a vertical temperature gradient which causes con- 
vection arises in moving media; in this case it is important to 
consider the interaction of convection with the shear velocity 
fields. 

Convection is not a unique process which performs the transport 
of heat in stars or planets. Frequently, along with convection, 
radiation involves the transport of a section of the heat flow, by 
a certain method of changing convection. 

Large density gradients in astrophysical convective zones 
Cthus, on the Sun the density changes by 6 orders of magnitude in 
the convective shell, and on Jupiter — by approximately 3-4 orders 
of magnitude) make it impossible to use a very simplified study of 
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convection with the Boussinesq approximation for incompressible 
liquids (see, for example, 111). 

In order to study the global convective phenomena under 
astrophysical conditions, it is necessary to consider the spheri- 
cal geometry of the layers in which it occurs. Finally, we must 
note one of the very important factors which complicate the cal- 
culation of the convective zones in astro- and geophysics: con- 

vection is always turbulent In them due to the enormous linear 
dimensions of these zones. 

Examples of special effects which complicate the convective 
motion In astrophysical objects may be given later on, but it is /5 
sufficient to illustrate the difficulties arising in a study of 
convection under such conditions. A more detailed discussion of 
the specific nature of astrophysical convection in stars and the 
Sun may be found in 12, 33. 

In this report, we discuss the influence of only individual 
factors upon convection, but in a more rigorous mathematical form- 
ulation. This influences the convection of spherical geometry. 
Secondly, it influences rotation and partially shear flow and 
the nonuniform distribution of temperature at the layer boundaries. 

A discussion of the influence of certain other effects may be 
found in the report by Spiegel [2] and the studies 13-83. 


Section 1. Development of convective instability in spherical 
layer at rest . 

Let us begin by examining the problem of convective stability 
of a viscous, incompressible liquid in a spherical layer of the 
thickness h, located between two concentric, spherical surfaces 
with the rad:'.i r^ ■ r Q h and r ■ (r 0 + l)h. A spherically symmet- 
ric radial gravitation field acts upon the liquid (below, we 
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shall use the spherical system ^ j V, tp ). The tempera- 
ture constants at the boundaries of the layers T^ and T g and 
the heat sources of constant Intensity q, which are distributed 
uniformly In the liquid, create a radial temperature gradient 
7T- When the gradient is small, the liquid is at rest, and 
the temperature distribution is determined by the heat conduc- 
tivity 

f77 ' = ~ 2F (fc + zF*) = 

' ' CL.1) 

re -T e - §ki Zr '+ * 

heat conductivity coefficient assumed to be constant. 

When VT is large enough, the heat conductivity mode becomes 
unstable, and convection arises. 

The linear problem of the development of convection in a 
spherical layer at rest was first examined in the studies of 
Chandrasekhar [4] . The dimensionless equations for small equili- 
brium perturbations in the Boussinesq approximation have the 
form: 

•ff ~ ~ Vp * Au. + RclF f(r)$ y 

Pr * u. yr * a 6 cLitru=0. fl 2) 

Here U,p and & — velocity, pressure and temperature of the 
perturbed motion, — temperature gradient determined by the 
heat conductivity; the length, time, veloc- 
ity, pressure and temperature are used as characteristic scales 
and are represented by 


U 


The following notation is introduced: 

Pr*W- - Prandtl number, 

the Rayleigh number; 

a^v and X — coefficients of volumetric expansion, viscosity 

and heat conductivity 

^ g 0 — density, temperature gradient, and accelera- 
tion of gravity at the external layer boundary, 

*(r) - *(r)/g. 

It is apparent that at the layer boundaries the following con- 
ditions are satisfied: 


at 


and 


rs r.+i 


( 1 . 3 ) 


If the boundaries are solid, then it is necessary that the remain- 
ing velocity components also vanish: 


Uj * Uy « Q. 


( 1 . 4 > 


The corresponding viscous tangential stresses must vanish at the 
free boundaries : 


m 0* 


( 1 . 5 ) 


Let us give the exponential dependence of the perturbations 

on time: U* f p , 9 ~~ erp St , 


a 


In order to exclude the continuity equation, it is convenient 
to use the general representation of a solenoidal velocity field 
in the form of the sum of the poloidal and toroidal components: 

u. * xot *sr + urtvr. 

( 1 . 6 ) 
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Applying the operations \ f*%ot and Fz,0t* 
to the first equation (1.2), we obtain: 

L*(A -6)w *0, 

L '(A-G)A s - Rol i(r) L *6 - 0. ° 

PreQ +u 7 T -ab =q. 

Here we use L? to designate the operator 

(/<•«*$). 


The boundary conditions in terms of the scalars of the 
poloidal § and toroidal W of the fields may be written in the 
form: 



on the solid boundary, 


9=S - =$?f~r’-0 


on the free boundary 


( 1 . 8 ) 


The spherical symmetry of the problem makes it possible to 
expand the solution in spherical harmonics 

t fs m ' 

( r\ * — associated Legendre functions). We set 

Mr. 3 v) -wfi) 3 ?(#, (f ) - etc. 


Then we obtain the following equation for radial perturbation ampli- 
tudes : 

6 


Cl. 9) 


fa-rjW --Q. ' 

$ e (SD t - s)$ - RcLt(r)9 *o , 

(\ -Prf>)9 + &^6(r)£ *0. 


(1 


where the differential operator D 0 is determined as 

© 





and the indices 1 and m are omitted here and below. 

Thus, the problem of the convective stability of equilibrium 
is reduced to the problem of the eigenvalues for the Rayleigh 
number Ra. 


It is known from (1.9) that W ■ 0, i.e», the convective flow 
arising at the stability limit is poloidai. We should note the 
following singularity of the problem being considered: The 

azimuthal wave number m is not Included in the boundary value 
problem (1.8) - (1.9). For this reason, the critical Rayleigh 
numbers will depend only on the wave number 1, but not on m, i.e., 
the following degeneracy holds: zt*i perturbations correspond 

to one and the same critical Rayleigh number with a different 

angular dependence 01(f) e ±tm * (msZ) . A situation 

arises which is similar to that which occurs in the plane layer. 

An infinite set of solutions corresponds to the same Rayleigh 
number, since Ra depends only on the square of the modulus of 
the horizontal wave vector , and the relationship 

between a x and a y is arbitrary fl,4]. 
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When the following condition is satisfied 


cL r X(r) a. 10) a 

aF*iro . 

the problem (1.8) - (1.9) is self-conjugate, and consequently 
the principle of the stability displacement §" • 0 is valid 
and the convection which arises at the stability limit is 
stationary. Numerical calculations later performed by Durney 
19] and Young 110] showed that the principle of the stability dis- 
placement is apparently valid not only under the condition (1.10), 
but also in the more general case of convection development in a 
quiet spherical layer. In the particular case \ /f r \ fLirto*** (r 
the system (1.9) assumes the form: ' ' ™. " » # ** 


$ - fto. (r. * /)" 9 *Q, 

=o. 


We shall discuss briefly the method of solving the system 
(1.11) with the boundary conditions (1.8), since it is character- 
istic for solving the more general problems of liquid motion in 
spherical layers. In many cases it is convenient to use the Galerkin 
method expanding the radial amplitudes in series with respect to 
the Chandrasekhar functions 



which represents a combination of Bessel 
order and which satisfy the conditions 


Etfaln*ojm0 


functions ofL 

Q*M*U 


jctions, o£ half-integral 
and 


where are the roots of the equation 


Ce (#i( r »+0) S Q* 


( 1 . 12 ) 


It is known that equation (1.12) has an infinite set of 
real simple roots, and the functions UPs) form a complete 
orthogonal system of functions with the integral propertv 





Perturbations of temperature and the toroidal component of the veloc- 
ity field, when it differs from zero, may be represented in the form 
of a series with respect to the Chandrasekhar functions, and the 
pololdal component of the velocity field in the following form: 


S Ct./i (*i r )A£r* ts/r '+ 

► • 




where the constants ^ t " t and are selected so as to 

satisfy the four boundary conditions applied to the function S. 

Then the problem may be solved according to the well-known procedure. 


The critical Rayleigh numbers for four types of boundary con- 
ditions are found in l^Jand given in the table, together with the 
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CRITICAL RAYLEIGH NUMBERS AND WAVE NUMBERS IN LAYERS OP DIFFERING THICKNESS 

AND FOR DIFFERENT BOUNDARY CONDITIONS 
















































values of the critical wave numbers. It may be seen from the table 
that layers with more rigid boundaries are more stable , and the 
convective cells in them have a smaller size than in layers with 
fr^fi boundaries.. A decrease in the layer thickness 

h/r.h * r~ leads (see the table) to a decrease in 
the dimensions of the convective cells, i.e., to an increase in the 
critical wave number 1. 


For a study of convection in the atmospheres of the planets 
and in convective shells of stars like the Sun, of greatest inter- 
est are the rules for convective flows in thin layers 0.25 

with free boundaries or with a rigid lower and free upper boundary. 


Within the limit of very thin layers at & d, ) 

the system (1.11) may be reduced to a system of equations for con- /12 
vection in an infinite horizontal layer [1,4], As is known, in 
this case for a layer with free boundaries, the critical RavleiKh 

£ A 2 HT m 

number Ra Q * the horizontal wave number * £ y/ 

~ J? 7 

and the radial amplitudes do not depend on the wave number 1 : 


£>( r 'sin T(r- r c ), Q(r) - % tin ’Ji (r- r c ). 


Section 2. Influence of rotation, shear, and latitudinal tempera- 
ture gradient upon convection in a slowly rotating 
layer 

For problems of astrophysics, there has been great interest in 
the development of convection in moving rotating media with non- 
uniform temperature distribution at the boundary. This formulation 
of the problem occurs, for example, when modeling the global hydro- 
dynamic movements in atmospheres of large planets, particularly 
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Jupiter 111). 


Actually, the existence of an internal heat source for Jupiter 
C12] must lead to convective instability of the lower part of its 
atmosphere and to establishment of a turbulent convection mode, 
since the Rayleigh numbers are enormous. 

However, the data from observations of the visible surface of 
Jupiter Cthe ordered band-like structure) point to the highly organ- 
ized nature of large-scale movements in the atmosphere of the planet 
Csee, for example, [13] ) . 

It must be assumed that large-scale convection having the form 
of axisymmetric rolls develops on the background of comparatively 
small-scale turbulent convection, which may be taken into considera- 
tion by the introduction of the effective viscosity coefficients. 

If the atmosphere of the planet is not comparatively dense, 
then the absorption of solar radiation plays an important role, 
which may cause both shear zonal flows in the atmosphere and heating 
which is not uniform with respect to latitude 111) . 

The study I11J examined convective stability of shear flow of 
a viscous liquid in a slowly rotating spherical layer with a lati- 
tudinal temperature gradient on the external boundary. For simplic- 
ity, it was assumed that the flow arises due to the difference in 
the angular rotational velocities and n(l + e 5 of the internal 
and external boundaries of the layer and the latitudinal temperature 
gradient. It is assumed that the temperature is constant on the In- 
ternal boundary. The latitudinal temperature gradient which simu- 
lates the heating of the planetary atmosphere by solar radiation 
is assumed to be symmetric with respect to the equator, i.e., the 
temperature on the external boundary may be represented in the form 
of a series in terms of even Legendre polynomials 


12 







T t (9) = Z. 

rt 





It is apparent that the temperature T g which is average 
the external boundary equals T^° . Let us examine the case 


on 




Due to two factors (differences of angular velocities of sphere 
rotation and latitudinal gradient), a stationary flow arises which 
may be described by the following dimensionless equations in a sy- 
stem of coordinates rotating at the angular velocity fi. 


(V v)v= -VP+AV + Pr'Ra i(r)r % 2 

' (2.1) 

Pr V VT = AT +CJ , ; cU<S V*Q 


with the boundary conditions 


r«r« : V sT -Q } 

r-r 0 +±; V= £ (#, t T- - /- 


( 2 . 2 ) 



I 

| 


! 

! 


1 


Here : 

T/.P and T — velocity, pressure and flow temperature. 

Re = — Reynolds number,^ a Tr^/T^^ 2* • 

unit vector along the axis of rotation. When deriving the equa- 
tions, it was assumed that the relationship of centrifugal force 
to the gravitational force is small, and VP includes all the 
potential forces. 


/14 


The stationary flow (2.1) - (.2.2) may be unstable with re 
spect to small perturbations at large e Re numbers — shear in 
stability, and at large Rayleigh numbers Ra — convective 
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instability. In the general case in the space of the parameters 

(e.fe , Ra, X ) we may find the surface of the neutral stability 
and determine the region of motion instability. 


The system of equations for small perturbations may be great- 
ly complicated as compared with the case of a fixed layer (1.7): 

-2 hQs - rvut [(a 7)V-f^7)o]-Q, 
[t i (d.- c ') + 2 fe^]^s +2J?eQ-w - Satyr) L*§* 
+?td t luV)U’+(Vr)u] = 0, 

Pr 6.8 * PrVVB +uVT -A8 * 0, «- 3 > 

where 

cc-D r‘- 

The boundary conditions remain the same (1.8). The dependence 
of the Coriolis force components on the angle in a spherical 
rotating layer does not allow separation of the variables in the 
system (2.3). However, if we are interested in convective stabil- 
ity of the flow in a slowly rotating layer for a small latitudinal 
temperature gradient at the external boundary, then this difficulty /15 
may be avoided by solving the system (2.3) with the boundary condi- 
tions (1.8) by the method of regular perturbations. 


The basic flow becomes unstable at a certain critical value 
of the Rayleigh number, which in the general case must depend on 
other dimensionless parameters: ^ 

If the characteristics of the basic motion and the perturba- 
tions are represented in the form of series with whole positive 




powers of Re and X [ill 




*0 


then the solution of the boundary value problem (2.3), (1.8) re- 
duces to solving a sequence of systems of equations with boundary 
conditions of the type (1.8). The system of equations for a zero 
approximation coincides with the system (1.7) for a quiet layer. 
Similarly, (2.U) may be represented in the form of the series 
RQ C and @ j , since the principle of stability displacement in 
rotating liquids is not valid in the general case and S CT 
differs from zero. We should note that the rotation, shear and 
latitudinal gradient, as will be seen below, remove the degeneracy, 
since the critical number Ra c depends on the azimuthal wave num- 
ber m and its minimization makes it possible to establish the form 
of the perturbations precisely which arise at the stability limit. 


It is apparent that the solution depends on the relationship 
between the parameters Re and A . Therefore, we set Re- X • 
Then in different ranges of changes in the parameter a , other 
physical factors will have the basic influence upon flow stabiliza- 
tion and the convection which is produced. 


In the case of a thin layer r 0 >>i and the free boundaries, 
the solution of the problem is obtained in analytical form for ar- 
bitrary values of a. The limited nature of the viscous stresses in 
a thin layer requires that £ , where _ /1 6 

It was shown in [11] that in the following intervals in which 
« changes 
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1. 0 *« < '/ z Gb ReS to , /?a c Bf?a,4.*e'*o*o, 

2 . ot * '/z , 6 "* Re6,i t ROf’xfKU, 

3 . 0 <*'/z ,6 * Re. 

qualitatively different solutions are obtained; Rcu — critical 
Rayleigh number for a spherical layer at rest. In the expansion 
of and M , only the main terms of the series are written 
which, in the given interval of change for a , take into account the 
contribution of the basic factor. 


At 0 4 Ot <* 1/2 the stability limit, the form and the phase 
propagation velocity of the critical perturbations w are determined 
by the rotation and the shear. Outside of the region of values 
for the parameters £^ and Pr, lying between the curves £ andc 0 

(Fig. 1, the convective flow has the form of axisymmetric stationary 

v 4 

rolls Cm ■ Cj). Between the curves c 0 and C 0 the convection is 
three-dimensional, has the form of banana-like cells with an azi- 
muthal wave number fy\~ £. , and is propagated in a wave-like manner 

in the y direction wltn the phase velocity a e, Re/zti. 

determined by the shear parameter . On the curves £ ?(Pr) 

the convection form is indeterminate, and it is necessary to con- 
sider the following expansion terms with respect to Re in order to 
eliminate the degeneracy. 


At £.>0 the critical Rayleigh number Rq c is always 
greater than Ra t , i.e., the shear flow in the rotation direction 
always stabilizes it independent of the form of the convection. A 
very large negative flow shear £ c 4~flo destabil- 

izes the convection, although the form of the convection does not 
change. The dependence of RQ t0 *R Q c-to o on the Rossby num- 

ber $ 9 is given in Fig. 2. 


The value £ m 0 corresponds to convection in a slowly rotat- 
ing but solid layer and was studied in detail in IlM], In this /17 

case, as may be seen from Fig. 1, convection is established in the 


16 


form of banana-like cells, but its velocity 

is always directed toward the Inverse rotation of the layer, and 
is much less than w Q , since it is proportional to r Q . We 
should note that in thin layers (greater than 1 ) the axisymmetric 
stationary convection Cm ■ Ql has maximum amplitudes in the polar 
regions, and on the other hand the banana-like convection Cm ■ 1„) 
is concentrated close to the equator. 

Allowance for the nonlinearity in this problem by studying 
the effects of second order terms of flow amplitude A shows 
that the redistribution of the angular convection moment leads 
in the second approximation to differential rotation of the layer. 
Actually, for the -components of the velocity field averaged 
over y and r, the following analytical expression was obtained 
in 114] 




*> = 

9t r 


A z Re 


£ /n & 


(P e ) W 


C2.51 


for the case of a thin layer and free boundaries. Numerical calcu- 
lations carried out for other boundary conditions show that the 
case of free boundaries to a certain degree is exclusive, since the 
value of the rotation greatly decreases for more general boundary 
conditions. However, in all cases the differential rotation is 
characterized by the equatorial acceleration of the layer. 


This model for the interaction of convection with rotation may 
he used to explain the observational data regarding the Sun. If 
the theoretically examined convective cells are identified with 
gigantic cells on the Sun, and the action of the convective flows 
of smaller scales is taken into account by the introduction of the 
effective viscosity, then the amplitude of the differential rotation 
in the case of free upper and solid lower boundaries according to /18 
Busse estimates [14] coincides with the observed equatorial accelera- 
tion of the Sun. 
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At C< - 1/2 , the total effect of shear, rotation and the 
latitudinal temperature gradient determines the parameters of 
critical convection. In Pig. 3, in the space of the parameters 
£ o,Pr. A /fie* , the surface JT given by the equation 

C ( Pr)£* - BC. -4+ 4 A /Re * *0 

* * 

at “ 8,986, Q m 2,49 Z, C(Pr)m 0 t SXPr t *O,2B9Pr * 0.87 

delineates, together with the planes A/R € • 0 and * 

the region within which the critical perturbations are realized in 
the form of three-dimensional banana-like cells with ptf • £ and 
outside of it in the form of the axisymmetric rolls with m ■ 0. 

The region of stabilization and destabilization of the basic flow 
is divided by the plane fj given by the equation 

R o * Be.- * /Re**0* 

Below this curve there is always stabilization of the basic 
flow; above it — destabilization. Axisymmetric convection is 
always stationary, and three-dimensional cells are propagated in 
a wave-like manner with the azimuthal phase velocity w Q . 

When Oi > y2 , the stability limit and the form of the 
critical perturbations are determined by the latitudinal tempera- 
ture gradient, and the decrement £ and the phase velocity w^ 
proportional to it — by rotation. Depending on the sign of X > 

the minimum critical number Ra„ will be achieved in the case of 

c 

axisymmetric convection, if /j > Q and three-dimensional cells 
if /(<•()• In Both cases > Ra c < Ra 0 » i*e., the nonuniform 
temperature distribution at the external boundary as compared 


18 


with the uniform distribution leads to destabilization of the basic 
mode for one and the same total supply of heat passing to the 
layer [15] . 


719 


Thus, the total action on convection of all three effects in 
thin spherical layers is determined by the values of the similarity 
parameters e Q , Pr and X/Re . In the case of slow rotation (small 
Re), critical convection is always three-dimensional and has a wave 
nature [14]. The shear and the latitudinal temperature gradient 
at X > o and in rotating layers coiiwr’.bute to the development of 
convective cells in the form of axisymraetric rolls. As applied to 
the astrophysical problem being considered, this means that allow- 
ance for solar radiation in a certain form may lead to a visible 
band-like structure of the planet surface. Actually, in adjacent 
rolls (two rolls form a convective cell), the gas in the meridional 
plane moves in opposite directions, forming ascending (warmer) and 
descending Ccooler] streams. If, as is assumed, we interpret the 
light bands (zones) as ascending streams in the atmosphere of 
Jupiter, and the dark bands (belts) as descending streams, then the 
picture obtained for axisyrametric convection qualitatively explains 
the visible band-like structure of the surface. 

If we assume that the band-like structure of the surface of 
Jupiter represents the occurrence of convective motions arising 
under the influence of an internal heat source and solar radiation, 
in the visible surface layer of the planet, and that this convection 
is. close to a critical value with respect to the effective transport 
coefficients, then on the basis of the results obtained we may de- 
termine the thickness of the convective layer and the effective 

g 

transport coefficients which are on the order of 150 km and 10 - 

7 2 

IQ cm /sec, respectively [11]. 
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Section 3. Influence of rotation on critical convection In 
the case of arbitrary Taylor numbers 

In the case of the rotation of a spherical layer with an art>i- 
trary velocity, the problem is greatly complicated, and complete /20 

separation of the variables cannot be achieved. It is possible to 
separate only the dependence in terms of longitude represented by 
the azimuthal wave number m. To study the influence of rotation 
on critical convection, it is necessary to solve a linear system of 
partial differential equations. 


Gilman ll61 obtained more complete results for this problem for 
a layer with the thickness 5 ■ 0.25 with free and infinitely heat- 
conducting boundaries. The calculations were performed for the 
Prandtl number Pr ■ 1. The system of linear equations in this 
Btudy was solved by the method of finite differences in the 
plane for given m from 0 to 24. A determination was made of the 
rate at which the perturbations increased and the frequency 
of oscillations of the most unstable symmetric and anti-syraraetric 
perturbations with respect to the equator for a given m. The Ra c 
numbers at the stability limit for each m and the structure of the 
perturbations were found by interpolation up to the zero rate of 
increase. The Taylor numbers m « *J2 V/y - 

changed fror 0 to 10 . The calculation results were the following: 


1) The rotation had the strongest stabilizing influence upon 
modes which were anti-symmetric with respect to the equator, and 
on modes with low azimuthal wave numbers m (Pig. 4 and 4a) ; 


2) With an increase in Ta, the wave number m c of the most un- 
stable perturbations increases, i.e., critical convection in a ro- 
tating spherical layer is always three-dimensional, and the dimen- 
sions of the cell in the azimuthal direction decrease with an in- 
crease in the rotation rate (Fig. 4); 
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3) With an increase in the rotation rate, the modes corre- 
sponding to different wave numbers m may be divided into two 
classes; a large part of the modes with high m have a maximum 

of amplitude on or close to the equator, and a small part with /21 

small m have a maximum at the pole (Fig. 5); 

4) The equatorial modes of both symmetries are distributed 

in the azimuthal direction toward a side which is Inverse to the 
rotation (w <0), for a alow rotation of the layer in accordance 
with the results of Busse II 1 * 1 and the side of rotation tw > 0) 
for large Ta (Fig. 6), Polar modes (with the exception of the 

stationary mode m ■ 0) are propagated at much slower rates; 


5) The critical convective modes for large Ta are character- 
ized by the asymptotic relations^^ ^ 7it ^ Tct?* CO^-T q,^ 

obtained earlier by Roberts [17 1 and Busse 132) for convection in 
rapidly rotating spheres; 

6) In rapidly rotating layers, symmetric equatorial modes 
assume the form of rolls twisted around an axis parallel to the 
axis of rotation. The antisymmetric modes also have the form of 
rolls, but the direction of motion in them has different signs in 
different hemispheres and the liquid is input through the equatorial 
plane parallel to the axis of rotation. With an increase in Ta, 
the rolls are pressed into the internal sphere, obeying the Taylor- 
Praudmen theorem. The polar modes turn in a circular vortex, press- 
ing closely against the poles; 

7) Determinations show that the radial heat flux has a maxi- 
mum at the equator for symmetrical equatorial modes, and at lower 
latitudes for antisymmetric modes. However, at large Ta, both 
fluxes are greatly suppressed close to the equator at the upper 
boundary of the layer. The symmetrical equatorial modes transfer 
the heat flux to the equator, and the antisymmetric ones — to the 
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pole at the lower latitudes and to the equator at the higher lati- 
tudes. The symmetric equatorial modes transfer the angular moment 
to the equator from the high latitudes for al values of the Ta 
numbers ; 

8) For layers of other thicknesses 5 * 0.11 and 0.67, the 
behavior of the critical convective perturbations is similar to 
that described. 

These results are confirmed by the studies of Roberts [ 17] /22 

and Busse l 32] on the asymptotic modes of critical convection in 

spheres for very rapid rotation (Ta " , Eckman number E ■ 

V i 

*j(p r <.<■ i » -adius of the sphere). The dynamics of 
the convection which can be realized physically must in this case 
obey the Taylor-Praudmen theorem, according to which slow, steady 
motion of a slightly viscous liquid cannot depend on the coordinate 
along the rotation axis. Since such motions cannot be realized 
either in spheres or in spherical layers, the convective movements 
predominate over the limitations imposed by the theorem due to the 
nonstationarity and increase in the influence of viscosity due to 
a decrease in the scale of motion in the azimuthal direction. 

At the stability limit, convection in the spheres arises close 
to a cylindrical surface located at a distance of approximately 
1/2 r c from the rotation axis, and has the form of thin cylindrical 
rolls parallel to the rotation axis and slowly moving in the azi- 
muthal direction at a velocity which depends on the Prandtl number. 

The horizontal scale of the rolls is approximately E (Pig. 7). 

It is clear from the form of the convection that the component 
of the gravitational force which is perpendicular to the rotation 
axis has the basic influence upon motion of this type. 

Apparently in thick spherical layers with convection » 

just as in a sphere, arises close to the cylindrical surface with 
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the radius r c /2, and in thin layers with hiK/z the convective 
rolls parallel to the rotation axis are tangent to the internal 
sphere, i.e. t they have a maximum possible vertical scale and, 
consequently, a minimum deviation from the Taylor-Praudmen theorem. 

Thus, at the stability limit in spherical layers both in the /23 
case of slow and rapid rotation, a non-axisymmetrlcal mode is es- 
tablished which is symmetrical with respect to the equator and 
with a high azimuthal wave number m. In the case of rapid rotation, 
the convection is concentrated in the pre-equatorlal band, and the 
phase velocity is directed toward the side of the basic rotation 
rate. 


Section Nonlinear convection in a quiet spherical layer 

The solutions of a linear system of equations which describe 
convection in a quiet spherical layer are, Just as in the plane 
case, degenerate. This means that the solutions with a different 
spatial dependence (in this case with different azimuthal wave num- 
bers m) correspond to one and the same eigenvalue — the critical 
Rayleigh number Ra . Their realization at the stability limit is 

v 

equally probable. This set of solutions 3s a reflection of the 
high degree of symmetry of the problem, but is this the physical 
result or is it due to disregarding Important physical effects? 

As is clear from the preceding, including even small additional 
effects (rotation, temperature nonuniformity, etc.) usually elimi- 
nates or decreases the degeneracy. One of the most important effects 
in the linear problem is the effect of the process nonlinearity. 
Therefore, it is important to establish whether allowance for non- 
linear terms makes it possible to determine, in addition to the am- 
plitude, the form of the finite-amplitude convection. 

Close to the stability limit (Ra - Ra Q ), the nonlinear terms 
may be regarded as a perturbation of the linear problem and the 
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asymptotic method may be used for expanding the solution in powers 
of the small parameter — the flow amplitude A. This method was 
used very successfully for studying the nonlinear convection 
motions Ra Q in the plane horizontal layer ( 1 83 . This method was 
used by Busse J 19 J to study the form of stationary finite-amplitude 
convective motions In a spherical layer. 

The nonlinear system of equations 


&L*s ~(iat(r)ft. 


F ?*{?*[ u * (V* U)]J, 


M * T'(r)r-'L*s * Pr u 79 


with allowance for the expansion in terms of amplitude 

t e-z, *% ,**•&*«' 


( 4 . 2 ) 


may be transformed into a sequence of linear nonhomogeneous systems 
of equations for ©^ and Ra^. The well-known system (1.9) is 
used as a system of the first order with respect to A. This system 
determines the critical convection in a spherical layer, the criti- 
cal wave number 1 and the Rayleigh number Ra 0 . 


The general solution of the system may be given in the form 


s, - S, (r) tr t ($ if) ; 6, * 6 } (rj *, ($. <f), 
ir € (<$, (o( m cosm<f+/i m sinm<?) 
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Here coefficients <*,. and B 

degeneracy of the system. 


remain arbitrary due to the 


The system of second order with respect to A may be written 
as follows 

A*L\ - RaJ(r)L*9 z « - L fflr 

-As t rVAs,]+ Rat(r)L*‘9 t ' 

A 9 x *i(t)rl\ = Pr[ V*( 7 *r$)]y 0 , . 


The solvability condition (4.3) and (4.4) has the form 

Ra., (is*L*&,) - <S *L '(Vfrr+ifo ns r AsSto$}><, 
- Pr Ra. (9 *[?(!+ r?)s, - r 'A S t ]r9 t >=Z7» 


where the brackets designate integration over the entire volume, 
and S + and 0 + — the solutions of a conjugate homogeneous sytem 

aY*s * - Ro. &(r)r~'L* 6 *- 0 , 

. , 0 . 6 ) 

A0 f t(r)L s* = 0< : 


It is convenient to represent the solution (4.6) in the form 

$ * = S *(r) if) , 9* -Q *(r) vf ( ’$,<?)> 

where (M) — spherical harmonics of order 1. Then (.4.5) 

may be written in the general form 


25 



Rci,<tr t *ir t >.-<Ui*tr 4 tr t >Mf ()•/?• 


(4.7) 


When the linear system corresponding to the nonlinear system 
(4.1) is self-conjugate, for any 1. For nonlinear 

systems of the more general form, the condition of the linear sy- 
stem self-conjugate property does not lead tr ^ s q 


We should note that the condition (4.7) is equivalent to the 
system it*t of nonlinear equations obtained due to the existence 
of lt*l independent harmonics of the order 1. 


The normalization condition v g supplements the last missing 
equation for determining 

z £ . Z i 

< Vk > = CXb +Z_ (°<m tp>m )=.£.. (4.8) /26 

Since the system (4.7) - (4.8) is satisfied by the value RQ t s O 

for any a and B for odd 1, the solvability condition does not 
mm 

eliminate the degeneracy of the solution for odd 1. It may be 
readily shown that all at 1 which are odd vanish. All even 

to* are usually positive, and we may thus conclude that in the 
case of odd 1 the subcritical instability ^ fa <Rai ) 
may not occur in the spherical layer. 


In the case of even 1, there is no general method for solving 
the system (4.7) - (4.8). Busse 119] examined several particular 
cases corresponding to 1 = 2, 4, 6 and very large even 1. A solu- 
tion of the system (4.7) - (4.8) for any 1 is the axisymmetrical 
solution 
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-0 , 
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with the corresponding value 
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In the case 1 = 2 the axisymmetric solution within an accuracy 
of the rotation transformation is unique and satisfies the condi- 
tion of solvability. 

In the general case of even 1 there are several solutions corre- 
sponding to different values of Ra^, i.e., the degeneracy is only 
partially removed by the solvability condition. A certain physical 
principle is therefore necessary for selecting the flow which can 
be realized in practice. However, usually for selecting the finite- 
amplitude flow, its stability is studied. However, this is a very 
complex problem, and Busse in this study confined himself to the 
"softest” principle, namely he assumed that close to the stability 
limit the finite-amplitude flows may be realized physically with 
a minimum value of the Rayleigh number, i.e., the flows correspond- /27 
ing to the point E in Pig. 7a. 

The minimum value of Ra is given by the formula 

where Raj is usually positive and reflects the change of the tem- 
perature field averaged horizontally due to convection, i.e., it 
changes very little with a change in 1. If we disregard these 
changes and the contribution of the terms of higher order, it is 
apparent that R a m ^ n will correspond to the largest value of Ra^. 

At 1 = 4 and 1=6, there are other stationary solutions in 
addition to the axisymmetric solution. Thus, in the case 1=4 the 
value of R is reached for solving with cubic symmetry (JPig. 8), 
and at 1 = 6 — for solving with dodecahedron symmetry (Pig. 9). 

For large even 1, only partial solutions are obtained. However, 
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the problem regarding the form of the most preferred mode remains 
open. 


If the corresponding values of Ra^fl) are not very small and J 

Ra Q (l) has no sharply expressed minimum for odd 1, then the realiza- 
tion of the convective flow with odd 1 is only slightly probable, | 

i 

and convective motions arise close to the stability limit in spheri- j 

cal layers as finite-amplitude subcritical instabilities correspond- 
ing to the even spherical harmonics. The form of these motions ap- 
parently has the greatest possible symmetry. f 

I 

For a very large supercriticality, expansion with respect to j 

amplitude is not applicable and it is necessary to solve the prob- 
lem numerically. Dumey [91 first attempted this, and calculated 
nonlinear convection in a spherical layer in the quas: linear ap- 
proximation. 

I 

The basic assumption used in this approximation is disregard- 
ing the nonlinear terms corresponding to the intrinsic interaction /2 8 
of the perturbations, i.e., only those nonlinear terms remained in 
the equations which describe the interaction of the average tempera- 
ture field with fluctuations of velocity and temperature [20] . 


As Herring [20] showed, the quasilinear approximation gives 
good quantitative results when calculating convective flows with 
large Prandtl numbers. In the nonlinear case, this approximation 
permits separation of the variables. The finite-amplitude solution 
may be represented in the form of one unique mode, when the solution 
is decomposed in series with respect to the spherical harmonics 


<fj . The system of equations of motion may be solved 
by the establishment method. Calculations for the layer with a 
thickness of 6 « 0.25 showed that convection at the stability limit 
and nonlinear motions in the case Rct.*5Q/ta a are always sta- 
tionary and correspond to the wave number 1 = 10, i.e., 1 equal to 
the critical value of the wave number for a given layer thickness. 


3 

i 
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Hie wave number realized In the case of convective motions may be 
determined by studying the stability of the solutions for glven^l 
with respect to small perturbations corresponding to and 

harmonics. For small supercriticality, the stable solution corre 
sponded to the solution maximizing the heat flow, in accordance 
with the Malkus principle [21]. 


Certain types of convection (for example, flow with cubic 
symmetry. Fig. 8) obtained by Busse analytically for small ampli- 
tudes were calculated by Young numerically for convective motions 
of arbitrary amplitude. 


Young [10] assumed a numerically complete nonlinear system 
of equations of hydrodynamics in the Boussinesq approximation for 
a spherical layer with free, infinitely heat-conducting boundaries. 

He used a rather original computational method: expansion in terms 

of spherical harmonics was used for angular variables, and the 
method of finite differences was applied with respect to the radius. 

The nonlinear terms were calculated using the modified transforms- /29 
tion method of Orszag [22], and the order of magnitude of the other 
derivatives with respect to r in the system of equations was reduced 
using the Green function for the operator and the one-dimension- 

al diffusion equation. 


The problem was solved by the establishment method for layers 
with the thickness 5 “2.33 and 0.67, the Prandtl number Pr * 1 and 
5 and the Rayleigh number Hat 5,5 Ra, Both axisymmetric and 
non-axisymmetrlc solutions were studied. 


Basic results: 


1) The axisymmetric solutions were not preferred for the com- 
plete three-dimensional equations of motion. There is an interval 
of values for the Rayleigi numbers ^ ^ 
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where there are simultaneously two steady state modes depending 
on the initial data: the axisymmetric and the non-axisymmetric. 

A similar result for the dependence of the solution on the Initial 
data was shown in the work by Schluter, et al. tl 8] lor a plane 
horizontal layer. Krishnamurti confirmed experimentally [23] that 
three-dimensional cells In a plane layer may exist in the region 
of values for the Ra numbers where the two-dimensional rolls are 
stable. 

2) At Ra > Ra(Pr,6) , the axisymmetric convection be- 
comes unstable with respect to three-dimensional perturbations. 

The ratio Ra/to. decreases with an increase in the layer thick- 
ness 5 and increases with an increase in the Prandtl number. 

Busse [24] obtained a qualitatively similar result pertaining to 
the instability of two-dimensional rolls with respect to three- 
dimensional perturbations in a plane layer. 

3) In contrast to the quasilinear approximation [93 , the 
dominating mode of the finite-amplitude solution does not always 
coincide with the most unstable mode at the stability limit. With /30 
an increase in the Ra number, in individual cases there is a sharp 
change in the horizontal flow structure (decrease in the dominating 
wave number 1), accompanied by a sharp change in the dependence of 

the heat flux on the Ra number (Pig. 10). This phenomenon has an 
analog in a plane layer and was observed experimentally in t23, 25, 

263 . 

4) The dimensionless heat flux determined by the Nusselt num- 
ber Nu Increases with an increase in Ra/Ra Q (Fig. IQ). However, 
the difference in the heat fluxes transmitted in axisymmetric and 
non-axisymmetric modes is small, on the order of 5 % (Fig. 11). 

5) _ For those values of the parameters when there are stable 
stationary modes, an oscillatory three-dimensional convection mode 
may be established depending on the initial data. The dependence 
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of the oscillation amplitude on time has almost a sinusoidal nature, 
but the oscillation periods of the heat flux and the poloidal and 
toroidal mode are different. 

It is important to note that toroidal modes in a quiet spheri- 
cal layer are always negligibly small, except for the case of the 
oscillatory mode. This fact and the absence of oscillatory modes 
in the quasilinear approximation indicates that the occurrence of 
oscillations is related to nonlinear terms which describe the in- 
trinsic interactions of the modes in the equations of motions. 

The results obtained by Busse 1 193 and Young HO] are insuffi- 
cient for formulating a general, although qualitative, representa- 
tion of nonlinear convection In a quiet spherical layer. On the 
other hand, the studies of Busse 119 J assume that close to the 
stability limit, just as in a plane layer when there is asymmetry 
between the upper and lower half of the layer, the three-dimensional 
cells are the preferred form of the finite-amplitude convective 
motions. In the opinion of Busse, the reason for this is the geometric 
asymmetry of the spherical layer. If we make an analogy with the 
plane layer later, then convection in the form of rolls must re- /31 
place the three-dimensional structure, when the amplitude of the 
motion becomes large as compared with the asymmetry. However, 
the numerical calculations of Young point to the opposite picture: 
the axisymmetric convection becomes unstable with respect to three- 
dimensional perturbations with an increase in the Rayleigh number 
Cat least for a Prandtl number which equals unity). 

However, we must keep in mind that the calculations of Young 
pertain to Rayleigh numbers R<X/Rq, 4 1.3 - 1.4, whereas the 

results of Busse H9J are limited to small supercriticality. It 
is possible that the case M(i) s o. which was not studied by Busse 
corresponds to a plane layer which is symmetric with respect to 
the middle, and convection in this case is excited very smoothly 
and has the form of axisymmetric rolls. 
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But this pertains to the area of analogies. Many problems 
are still not explained. Is there an axisymmetric solution with 
a unique stable solution for small supercriticality? If the 
answer is yes, then in the space of the parameters (Pr,6 ) where 
does the stability boundary lie of axisymmetric convection with 
respect to three-dimensional perturbations? What will occur in 
very thin layers 6 -*-0 (correct limiting transition to plane layer)? 

What is the influence of initial perturbations on the establishment 
of a certain form of convection motions? What is the influence of 
the Prandtl number on the transition from stationary convection to 
oscillatory convection? 

Section 5. Nonlinear convection: interaction with rotation. 

Results of numerical calculations. 

Calculations of the nonlinear convection in a rotating spherical 
layer were first performed by Durney [27j in the quasilinear approxi- 
mation. This approximation in the case of allowance for rotation 
makes it possible to greatly simplify the problem due to the fact 
that modes with different azimuthal wave numbers m in the equations /32 
do not interact directly. The case of azisymmetric convection 
m ■ Q was examined for a layer with the thickness $ * 0.25, the 
number Ra * 1500 and the Taylor number Ta < 500. The main influence 
of rotation on axisymmetric convection is reduced to its strong 
stabilization in the equatorial region, where pulsations of velocity, 
temperature, and the convective heat flux H*'/i UiB.iinS 
are suppressed (Pig. 12a). 

As the calculation showed, an increase in the angular rotation 
velocity requires considering a large number of harmonics for an 
adequate description of convection. Thus, in the case Ta » 500, 
close to the value Ta for which convection is completely suppressed, 
it is necessary to consider modes with wave numbers from the interval 
20 £ 6 £ 6 * 
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However, it is- known that the strongest stabilizing influence j 

of rotation is exerted on axisyrametric convection [28, 291. There- 
fore, it is important to examine the non-axisymmetric case m s* 0, 
which was done by Durney [30 ] for the number Ra ■ 1500 and the nura- f 

her Ta » 4. 

To establish the preferred form of convective motion, integra- 
tion was performed over time of a system of equations in which 
several modes were retained corresponding to different m. With 
time, all of the modes were damped except m - 10. In the equations, j 

modes were retained with poloidal wave numbers 1 ■ 8, 10, 12 for 
temperature and the poloidal field component and 1 - 9» 11 for the 
toroidal component. We should recall that 1 » 10 is the critical 
wave number for a spherical layer at rest with the thickness 
6 « 0.25 and the mode l c - 10 is the unique stable mode when cal- 
culating nonlinear convection in a spherical layer at rest in the 
quasilinear approximation at Ra * 1500. 

Thus, for small Taylor numbers, the convective flow arises in /33 
the form of three-dimensional cells with m = 1 , which coincides 
with the results of Busse for critical convection in thin layers 
Il4], The dependence of the solution on time has, in accordance 
with [l4j, a specific form expL(my>4$t), which is characteris- 
tic for gravitational-hygroscopic waves. The value of does 

not depend on r and Q , and in this case equals 3.755. 

The convective heat flux C Cl & } » averaged over Ip , does 
not depend on time and is only a function of t and $ . Figure j 

12b shows the dependence <u.9> on in the middle of the j 

layer. It may be seen that the heat flux is maximal at the equator, j 

i.e., non-axisymmetric convection is more effective at the equator, 
in contrast to axisymmetric convection. 

After the solutions are found in the quasilinear approximation, 
we may calculate the nonlinear terms corresponding to the intrinsic 
interaction of the modes and determine how the angular moment is 


.i 
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redistributed by convective motion. This procedure is valid if 
the nonlinear terms are small as compared with terms which are 
considered in the quasillnear approximation. Figure 13 shows 
the dimensionless part of the angular frequency depending on /$, 
at the external boundary as a function of the angle 


Thus, although the convective motions are non-axisymmetric, 
the intrinsic Interaction of non-axisymmetric modes in the subse- 
quent approximation excites the axisymmetric mode, leading to dif- 
ferential rotation of the layer and equatorial acceleration. 


This property of convective motion causes differential rota- 
tion at the external boundary and may be used to explain the 
equatorial acceleration of the Sun (30J. 
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The studies of Gilman [ 16 ] and Busse tl 1 !] of critical convec- /3 4 
tion in a rotating spherical layer and of Young [10] and Busse 
[19] on nonlinear convection in a quiet layer showed that both 
effects (rotation and ronlinearity ) operate in one direction, i.e., 
they cause non-axisymruetric convection. Durney came to the same 
result 130], when studying the interaction of nonlinear convection 
with rotation in a quasilinear approximation. 


Therefore, at first glance the result obtained in the study 
of Williams and Robinson r 313 is somewhat unexpected, where a 
numerical calculation was made of nonlinear convective motion in 
the approximation of a "semi-thin layer" for thefts < 9502a 9 RQ q 
number and the Taylor number 0 £ TCL $ 5 # IQ? and 

Pr ■ 7.1 (Jta Q — critical Rayleigh number in a quiet layer). The 
calculations employed the grid method, the equations were integrated 
over time for a layer with the thickness 5 ■ 0.15 with a free upper 
and rigid lower boundary, and both boundaries were assumed to be in- 
finitely conducting. 

The basic result is reduced to the fact that in this range of 
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Ra and Ta numbers, a stationary axisymmetrlc convective mode is 
established which is stable with respect to three-dimensional 
(azimuthal) perturbations. 

It may be assumed that this divergence in the results may be 
explained by the large value of the Prandtl number used in [31 ]• 
However, this is not the case and may be seen from calculating 
A-12 (see [31 1, table). Actually, this case Ra ■ 2376 and Ta ■ 
2.9*10^ pertains to convection at the stability limit (Nusselt 
number Nu - 1.001), and as is known, the Prandtl number is not 
included in linear equations in the steady state case. The deci- 
sive factor is the fact that the authors used the ’’semi-thin” 
layer approximation, in which they disregarded the components of 
the Coriolis force, proportional to Siftv and making the greatest 
contribution at the equator . (i.e . , the r-component and the part ( 

SLu^sinV — the azimuthal component of the Coriolis force). 

The authors explain the difference with the results of calculat- 
ing critical convection in a thin layer Il4] by the fact that the 
analysis of Busse applied only to slow rotation. 

The behavior of the axlsymmetric solution 131) differed great- 
ly from the axisymmetric flows of Oilman tl6) and Dumey 127). In 
contrast to the latter, the axisymmetric flow in the approximation 
of a ”semi-thin” layer was suppressed at the poles, and not at the 
equator with an increase in the angular velocity of rotation (Fig. 
14). 

As may be seen from the latitudinal velocity profiles (Fig. 

14), for different Ra and Ta, the zonal flow in the solution which 
is symmetric with respect to the equator has a maximum close to 
the equator at the upper layer boundary. There is a small de- 
crease in velocity at the equator itself. 

The axisymmetric structure of convective flow, the accelera- 
tion of the flow close to the equator, the alternating ascending 
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and descending flows in adjacent axisynanetric rolls make this model 
of convection very interesting for explaining the observed phenomena 
on the surface of Jupiter. 



The convection calculation performed by the authors, taking 
into account the planetary scales, showed that equatorial accelera- 
tion corresponding to the observational data may be obtained for 
a comparatively fine atmosphere of approximately 150 - 500 km. 
However, it must be assumed that there is strong anisotropy in the 
transport coefficients in the horizontal and vertical directions 




Section 6. Laboratory modeling of astrophyslcal convection. 


The difficulties arising when creating a gravitational field 
with spherical symmetry under laboratory conditions 'ead to the 
fact that the experimental studies on convection are limited to 
a plane layer. However, convection in a plane layer has the prop- /36 
erty of horizontal isotropy, which does not exist in rotating 
spherical layers with Intrinsic gravitation. Therefore, the anal- 
ogy with experiment is very poor. The greatest difference is due 
to the fact that the angles between the force of gravitation and 
the angular velocity vector in plane and spherical cases do not 
coincide. 


However, in the case of rapid rotation, as shown by Busse 
133 J» due to the specific form of convection in the form of columns 
orrallel to the axis of rotation, the basic force is the component 
of the force of gravity perpendicular to the angular velocity. 

This fact is the starting point for laboratory modeling of astro- 
physical convection performed in the division of Space Physics 
of Columbia University 134 ]. The force of gravitation was modeled 
in the experiments by centrifugal force arising during the rotation 
of a spherical layer around a vertical axis. Figure 15 shows a 
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diagram of the laboratory equipment. To compensate fcr the dif- 
ference in the signs of centrifugal and gravitational force under 
laboratory conditions, the direction of the temperature gradient 
was changed to the opposite direction. The spherical layer be- 
tween the external (transparent) spherical shell and the Internal 
sphere was filled with liquid, to which aluminum flakes were 
added to visualise the flow. The layer was placed in a trans- 
parent container, in which thermostatically controlled water cir- 
culated to maintain the temperature at the external boundary of 
the layer. To maintain a constant temperature of the internal 
sphere, the thermostatically controlled water was added along the 
axis of rot et ion at a lower temperature. 
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The preliminary results now published qualitatively describe 
certain convective modes which are produced. Por a high rotation /37 
rate and a corresponding AT , a very non-axisymmetric convective 
mode was established, having the form of vertical columns regularly 
distributed in space (Pig. 7a), With a further increase in the 
rotational velocity, the convective columns were filled with the 
external section of a spherical ring. Due to the difference in the 
phase propagation velocities, the location of the columns in space 
was irregular. However, their form was rigorously parallel to the 
rotation axis. 


Hie thin layers make it possible to implement flows with rela- 
tively low Taylor numbers without decreasing the rotational velocity 
below the limit when the gravitational force begins. In this case 
convection assumes the form of banana-like cells obtained theoretic- 
ally in slowly-rotating thin layers C 1^1 . 

However, there is another possibility ior modeling convection 
in spherical layers with central gravitation The basic idea is 
to produce a central field of forces by an eleccric field. It was 
tested on equipment with spherical geometry by Chandra and Smylle 
l 35l at the Columbia University in Canada. 
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In this case the central field of forces was produced by the 
difference in the potentials y on external and Internal spheres. 
Heating the Internal sphere and cooling the external one estab- 
lished the necessary temperature gradient. The calculation showed 
that the volumetric force directed along the radius of a spherical 
layer has the value 
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Where", V*r' - temperature coefficient of the dielectric con- 
stant. For a thin layer 
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In dielectrics of the silicon oil type at T 

(^) 0 ~ *1tO “C ' QH m .JO 
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In order that the force of the Earth attraction does not have 
a great Influence upon the convection produced, the following is 
necessary: 


y,/^ » £. 

If we set f (£ » 10, then for a field strength of •ft 

■ 100 kV/cro, the following condition is obtained for the geometry 
of the layer • 7-3 cm, i.e., the dimensions of the equlp- 
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ment a}?e very limited. The basic limitation i*i imposed by the 
strength of the field r/h . It must not exceed a certain value 
of *100 kV/cm, at which breakdown may occur in the liquid. 


It is apparent that the use of the electric field for model- 
ing convection in a spherical layer was implemented most effec- 
tively by Busse and Carrigan, but the equipment was more complex 
to use. 


CONCLUSION 

The problem of convection in a spherical layer without con- 
sidering any complicating factors (rotation, shear, etc.) is de- 
termined by three-dimensionless parameters: the Rayleigh number 

Ra, the Prandtl number Pr, and the dimensionless thickness of the 
layer $ t i.e., one parameter of similarity more than in the prob- 
lem of Raleigh-Benar . 

Just as in the problem of convection in a plane layer, the 
linear problem is degenerate, but solutions correspond to 

the critical Rayleigh number in a spherical layer, and not an in- 
finite set, and ^ greatly depends on the layer thickness and, 
to a lesser extent, on the boundary conditions. 


Allowance for small supercriticality of the motion does not 
completely eliminate the degeneracy, although it greatly reduces 
the possible solutions H9] • Just as in the plane case, for one 
and the same values of the similarity parameters, there is nor- 
uniqueness of the solutions for nonlinear convection. However, in 
contrast to the plane layer, the horizontal dimensions of finite- 
amplitude flow are completely determinate, and the arbitrary selec- 
tion of a poloidal wave number is impossible. This is due to the 
closed nature of the volume. The nonuniqueness of the solution 
may be due to the form (azimuthal wave number) or nonstationary 
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nature of the solutions. 

The existence of such complicating factors as rotation, lati- 
tudinal temperature gradient at the boundary, shear flow leads 
to the elimination of degeneracy in the linear problem [111 . 

Shear and the latitudinal gradient with %>n lead to establish- 
ing axisymmetric convection, and rotation causes three-dimensional 
cells having an intrinsic phase propagation velocity in the azi- 
muthal direction. In thin layers, even in the case of slow rota- 
tion, convection is concentrated in the equatorial region. 

At large rotation velocities , critical convection 

in the polar regions is suppressed, and convective cells assume 
the form of thin rolls — columns parallel to the axis of rotation. 
The critical convective modes at Tcl-*0o are characterized by the 
asymptotic relations! 


I 


! 

I 

i 


" Rcl c ~Tq.™ m e ~Ta.' At ) ~7a* 

A study of the interaction between nonlinear convection and 
rotation has only begun. The studies of Durney [27, 30 J and /40 

Williams and Robinson [31J pertained to specific approximations 
which cannot give a complete representation of nonlinear convec- 
tive flows in spherically rotating layers. The difficulties con- 
sist of the fact that modes of nonlinear convection which can be j 

physically realized (stable) in this case are always three-dimen- * 

sional, which greatly complicates the use of numerical methods. j 

The studies performed on convection in spherically rotating 
layers show the promising use of the simple models for understand- 
ing and explaining the observational phenomena on planets and stars. 

These models help to explain the most characteristic aspects of 
the phenomenon, in the case of convection in rotating spherical 
layers. These models also provide a basis for understanding the 
mechanism of the phenomenon as a whole and the formulation of more 
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complex models which are closer to real conditions, whose numerical 
modeling is still very difficult [33*36]. 

Along with this, certain difficulties have arisen. As was 
shown in [14, 30, 3H , the interaction of rotation with convection 
leads to two basic effects: transport of the angular moment to the 

equator and stabilization of large-scale convection at the poles. 

A consequence of the second effect is the great difference in the 
heat fluxes and temperatures at the poles and the equator at the 
upper boundary of the convective zone; this difference was not ob- 
served either on the Sun 1373 nor on Jupiter 138]. 

In addition, convection in rapidly rotating spherical shells 
strives to assume a non-axisymmetric form in the form of cells ex- 
tended in the meridional direction, and this contradicts the ob- 
servations on Jupiter. 

Recent results of spectroscopic measurements on "Pioneer 10 /4l 

and 11" 138] and theoretical studies 133, 36] show that the convec- 
tive model of Jupiter must be more complex than was assumed previous- 
ly. 


The model of Jupiter must explain the absence of a difference 
in the heat fluxes at the equator and at the poles and the develop- 
ment of convective cells in the form of axisymmetric rolls parallel 
to the equator. Studies we performed earlier 111, 363 lead to the 
assumption that both of these phenomena may be related to the ab- 
sorption of solar energy close to the upper boundary of the con- 
vective zone. 

Actually, it is difficult to assume that the heat flux from 
the Sun, which comprises approximately half of the internal heat 
flux, does not have a great influence upon the hydrodynamics of 
large-scale phenomena in the atmosphere of Jupiter. Apparently, 
in the upper layers of the atmosphere where primarily solar radia- 
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tion is absorbed, the viscosity is small (the penetrating convec- 
tion is suppressed due to the stable temperature stratification) 
and due to a geostrophic balance zonal westerly flows arise. The 
action of these flows on the underlying layers may be represented 
as tangential stresses jp , which are continuously distributed 
at the boundary and which change with latitude. In the case of 
rapid rotation of the planet Re and for small Rossby num- 
bers £ , the influence of on the lower layers 

leads to the formation of a cylindrical shear layer parallel to 
the axis of rotation and tangent to the lower boundary of the con- 
vective zor ^t the equator. As was shown in I36j, the extension 
of this layer to the external boundary causes a sharp change in 
the angular velocity at the surface, and may thus explain the 
existence of an equatorial stream in the region of + 10° of the 
Jupiter latitude. 

In addition, a strong zonal flow at the upper boundary, 
together with the latitudinal temperature gradient IllJ , in the 
presence of rotation contribute to the development of convective /42 
cells in the form of axisymmetric rolls. The band-like structure 
of the surface of Jupiter at latitudes of + 10° - + 45° may be the 
result of large-scale convective cells of this type reaching the 
surface . 

The difference in the absorption of solar radiation at the 
equator and the poles must lead to a decrease in the convection 
intensity in the pre-equatorial zone and may balance the influence 
of rotation causing a decrease in the intensity of convection at 
the poles. Actually, both of these factors, rotation and absorp- 
tion of solar radiation, will suppress the convection selectively: 
rotation — at the poles; absorption — close to the equator. As 
a result, the heat flux produced by the planet at different lati- 
tudes may be approximately the same. 
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The three-dimensional structure of convection at high latitudes 


(38] must be explained by a decrease in the velocity of the zonal 
flow and the absorption of solar energy 111] . It is apparent 
that the model approximates the real one most closely if convec- 
tion is regarded in a compressed medium, in contrast to the 
Boussinesq approximation. 

The situation is more complex with the solar convective model. 
Although the model of the convective shell has good agreement with 
the observational data referring to large-scale motions in the 
photosphere (differential rotation, dimensions of gigantic cells, 
weak meridional circulation, etc.), there is a great contrast be- 
tween the observed and calculated heat fluxes on the surface. 

Rough estimates of the difference in the heat fluxes at the 
equator and the poles &P for a deep model of the solar convective 
zone give values of /P . There is no 

validity to the ideas that the equations and boundary conditions 
are not sufficiently realistic for predicting heat fluxes close to 
the upper boundary of the convective zone, since the mechanism of /43 
transfer from convective heat transport to radiant heat transport 
close to the boundaries was not considered. Additional concepts 
regarding the possible redistribution of heat flux within the con- 
vective zone are necessary. 
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Regions corresponding to different forms of convective 
flow in the plane of the parameters (£•# )• 

Dependence of Ragg on the Rossby number for differ- 
ent values of the Pr numbers. 

Regions corresponding to different forms of the convec- 
tive flows in the space of the parameters 

( e. 

The Rayleigh numbers corresponding to the onset of in- 
stability which is symmetric with respect to the equator 
as a function of the wave number m for different Ta. 

The crosses designate critical modes. 

Ratio between the critical Rayleigh numbers for anti- 
symmetric and symmetric modes as a function of the Ta 
number. 

Normalized amplitudes of the mean square velocity aver- 
aged over the radius for ra ■ 0 - 24 as a function of the 

latitude, at 72t - 0, HJ* , 3-10* and HJ*. 

Dimensionless frequencies of oscillations of symmetric 
modes at the stability limit for different Ta. The 
crosses designate critical modes. 

Qualitative drawing of convective flows in rapidly rotat- 
ing layers. 

Schematic dependence of flow amplitude on the Ra number. 
Form of convection with cubic symmetry in the case 1 ■ 4. 
The solid lines — lines of constant radial velocity. The 
direction of motion is different in the cross hatched and 
non-cross hatched regions. 

Form of convection with dodecahedron symmetry in the case 
of 1 ■ 6. 

Nusselt number Nu as a function of in the case of 

axisymmetric convection at Pr ■ 5 for layers of differing 
thickness . 

Nusselt number as a function of Ra/Ra Q for axisymmetric 
(I) and non-axisymmetric (II) modes. 
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Pig. 

Pig. 

Pig. 


Pig. 


12. Convective heat fluxes for the number Ra ■ 1500 in the 
middle of a layer with the thickness £ ■ 0.25 as a 
function of the angle 0 

a) axisymmetric mode M ■ 0, Ta ■ 500 j 

b) non-axisyrametric mode m ■ 10, Ta ■ M. 

13. Dimensionless part of angular velocity depending on % 
at the external boundary of the layer. 

14. Latitudinal profiles of normalized zonal and radial 

w velocity and temperature in the middle of the layer 

for different values of Ra and Ta. 
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15. Diagram of equipment for modeling the astrophysical con- 
vection. 
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